arXiv:astro-ph/0503682v2 25 May 2005 


Mon. Not. R. Astron. Soc. 000.Q]-?? (2002) Printed 2 February 2008 (MN FTgX style file v2.2) 


The cooling function of HD molecule revisited 


Anton Lipovka 1 *, Ramona Nunez-Lopez 1,2 , and Vladimir Avila-Reese 2 

1 Centro de Investigation en Ft'sica, UNISON, Rosales y Blvd. Transversal, Col. Centro, Edif. 3-1, Hermosillo, Sonora, Mexico, 83000 
2 Instituto de Astronomta, U.N.A.M., A.P. 70-264, 04510, Mexico, D.F., Mexico 


2 February 2008 


ABSTRACT 

We report new calculations of the cooling rate of primordial gas by the HD molecule, tak¬ 
ing into account its ro-vibrational structure. The HD cooling function is calculated includ¬ 
ing radiative and collisional transitions for ./ ;C 8 rotational levels, and for the vibra¬ 
tional levels v = 0,1,2 and 3. The ro-vibrational level population is calculated from the 
balance equation assuming steady state. The cooling function is evaluated in the ranges 
of the kinetic temperatures, 7';., from 10 2 to 2 x 10 4 K and the number densities, rt//, 
from 1 to 10 8 cm -3 . We find that the inclusion of collisional ro-vibrational transitions in¬ 
creases significantly the HD cooling efficiency, in particular for high densities and tem¬ 
peratures. For tin 10 5 and 7), ~ 10 4 K the cooling function becomes more than an 
order of magnitude higher than previously reported. We give also the HD cooling rate in 
the presence of the cosmic microwave radiation field for radiation temperatures of 30, 85 
and 276 K (redshifts of 10, 30 and 100). The tabulated cooling functions are available at 
http://www.cifus.uson.mx/Personal_Pages/anton/DATA/HD_cooling/HD_cool.html We dis¬ 
cuss the relevance to explore the effects of including our results into models and simulations 
of galaxy formation, especially in the regime when gas cools down from temperatures above 
~ 3000K. 

Key words: cosmology: first stars — galaxies: formation — molecular processes 


1 INTRODUCTION 

The cooling and thermal balance of the primordial, zero-metallicity 
gas by molecular hydrogen (H 2 ) and deuterated hydrogen (HD) are 
key ingredients in the formation process of the first baryonic objects 
in the universe, such as dwarf-sized galaxies and Pop-Ill stars (see 
for recent reviews on this topic Barkana & Loeb 2001: Ciardi & 
Ferrara 2004; Bromm & Larson 2004). It is well known that H 2 and 
HD molecules form in the universe after recombination (e.g., Lepp 
& Shull 1984; Puy et al. 1993; Palla, Galli, & Silk 1995; Galli & 
Palla 1998; Stancil, Lepp, & Dalgarno 1998). The formation of H 2 
and HD molecules in primordial gas is also possible in a postshock 
flow (Mac Low & Shull 1986; Shapiro & Kang 1987; Uehara & 
Inutsuka 2000), where the gas density and temperature, as well as 
the molecular fractions, are different from those of the expanding 
homogeneous gas. 

The H 2 molecule is considered typically as the main coolant 
in the primordial medium (e.g., Tegmark et al. 1997; Abel, Bryan 
& Norman 2000; Bromm, Coppi & Larson 2002). Thus, rather 
detailed calculations were carried out for it (see for the most re¬ 
cent results and references therein Le Bourlot, Pineau des Forets & 
Flower 1999; Flower et al. 2000; Shaw et al. 2005). Nevertheless, 
in the case of dense primordial gas, the role of the HD molecule 
in the thermal and cooling process may become comparable to or 
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more important than the one of H 2 . Therefore, detailed calculations 
of the HD cooling function (the rate of cooling per HD molecule, 
hereafter CF) are needed. 

The HD molecule is an efficient coolant in the primordial 
medium because it has a permanent electric dipole moment that al¬ 
lows high probabilities for the radiative rotational transitions from 
the rotational levels pumped by collisions with H and He. Besides 
the HD molecule was efficiently produced in the reaction of chem¬ 
ical fractionation D + + H 2 —> H + + HD, which leads to a 
strong enhancement of the initial abundance of HD, up to roughly 
[HD]/[H 2 ] = 1CT 2 - 10“ 3 (e.g. Puy et al. 1993; Galli & Palla 
1998; Flower 2000). The HD molecule can cool the primordial gas 
to temperatures below ~ 300 K because it has allowed dipole rota¬ 
tional transitions characterized by energies two times smaller than 
those of the quadrupole transitions of H 2 . In that regard higher tem¬ 
peratures (> 3000K), cooling by HD can be again as important as 
that by H 2 . This question will be discussed in the present paper. 

The HD CF was calculated for the first time by Dalgarno & 
McCray (1972) for low temperatures and Boltzmann distributed 
rotational level populations. Varshalovich & Khersonskii (1976) 
improved the calculations by taking into account the departure 
from local thermodynamic equilibrium. More recently, Flower et 
al. (2000, and see also more references therein) have calculated the 
HD CF for low and intermedium kinetic temperatures and a large 
range of densities. These authors considered the rotational transi¬ 
tions within the vibrational ground state (v = 0) of HD. Note that 
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in the case of higher kinetic temperatures (Tk > 2 x 10 3 K), the 
populations of the vibrational levels become significant, and cool¬ 
ing by ro-vibrational transitions of HD may be important due to (i) 
the relatively high energy of vibrational quanta, and (ii) the rela¬ 
tively high probabilities of the ro-vibrational radiative transitions. 
Therefore, the inclusion of the the vibrational structure of the HD 
molecule in the calculation of its CF might be important. 

The aim of this paper is to recalculate the HD CF for a wide 
range of kinetic temperatures and densities by taking into account 
ro-vibrational transitions within ground (v = 0) and the first three 
excited (v = 1, 2, 3) vibrational levels. Further exploration of the 
influence of the HD CF calculated here on the thermal evolution 
of the primordial medium will be presented elsewhere. In §2 we 
give details of our calculations. The revised HD CF is given and 
compared with that from Flower et al. (2000) in §3. Finally, in §4 
we discuss our results and speculate about possible implications. 


2 THE MODEL 

It is well known that in the low density limit (uh < 10 3 cm -3 ), 
the relative populations of the levels of HD molecule depart sig¬ 
nificantly from the Boltzmann distribution. Therefore, to calculate 
the HD CF in the general case, we need the correct values of the 
ro-vibrational level populations outside the thermodynamical equi¬ 
librium. This can be calculated by using the equations of detailed 
population level balance: 

n v j + W vJ^VJ-) = 

V 1 J' 

V 1 J' 

^2,n vJ = 1 , ( 1 ) 

vJ 

where n v j is the population of the ro-vibrational level vJ, 
and are the probabilities of the radiative and 

collisional transitions, respectively. We assume a steady state for 
the population levels. This assumption has been used and discussed 
to be valid in our context by many other authors (e.g., Flower et al. 
2000; Flower & Pineau des Forets 2001; Le Petit, Roueff & Le 
Bourlot 2002). 

The population of the ro-vibrational levels of the HD 
molecule was calculated in the wide range of number densities nn 
from 1 to 10 s cm -3 , and for the kinetic temperature Tk from 10 2 to 
2 x 10 4 K. The ro-vibrational radiative transition probabilities for 
the HD molecule were taken from Abgrall et al. (1982). Regard¬ 
ing the collisional transition probabilities, it was shown in Flower 
et al. (2000) that the HD CF is insensitive to the H/H 2 density ra¬ 
tio. Therefore, one may take into account the excitation and de¬ 
excitation of HD only by collisions with the H atoms. The proba¬ 
bilities of the collisionally excited pure rotational transitions of HD 
(up to the third vibratioanally excited state, 0 ^ v ^ 4, and up to 
the 8-th rotational level 0 ^ J ^ 8) were taken from Roueff & 
Flower (1999) and Roueff & Zeippen (1999). We consider only the 
dipolar transitions (A J = ±1). Owing to the permanent dipolar 
moment of the HD molecule, these transitions dominate collisional 
population transfer (Flower & Roueff 1999). 

Flower & Roueff (1999) computed probability coefficients 
for several ro-vibrational transitions in HD induced by collisions 
with atomic and molecular hydrogen in a range of temperatures 
from 100 to 2000 K. In that paper, some values of the transition 
probabilities for the vibrational relaxation v = 1 — * 0 are 


presented and compared with the corresponding coefficients 
for H 2 . The rate coefficients for vibrational relaxation of HD 
and H 2 in collisions with H are similar in magnitude for a 
given temperature. Besides, the probabilities for both molecules 
change approximately in the same way with temperature. On 
the basis of these results and in the absence of extensive cal¬ 
culations for the HD+H collisional ro-vibrational transition 
probabilities in the range of temperatures we want to explore, 
we will use the corresponding probabilities of the H 2 +H (Tine, 
Lepp & Dalgarno 1998) given in electronic form on the website 
http://www.physics.unlv.edu/astrophysics/h2h2rates/index.html 
(Lepp, Tine & Dalgamo 1997, in preparation). 

The total CF per unit of volume is defined, for example for the 
HD molecule, as: 

A hd = uhdWhd, (2) 

where uhd is the HD number density, and Whd is the HD CF in 
unities of erg s _1 per HD molecule: 

w hd= ^ {n v jW R }^ v , J ,-n v ,j l W R \,^ vJ )hv v j^ v :j l .{'i) 

vJv' J' 

Here, W R j= A Ri - + B R> - u are the probabilities of the ra¬ 
diative transitions with the emission of the photon hu v j_ v i ji , and 
W„i ji _ v j — B Rl u are the probabilities of radiative transition 
with the absorption of the corresponding field photon. The sym¬ 
bols A and B are for the corresponding Einstein coefficients and 
u is the radiation field. The populations of the ro-vibrational levels 
n v j are those calculated from the balance equation (1). This defi¬ 
nition of Whd is more general than the one used in Flower et al. 
(2000), because it takes into account the potential effects of a ra¬ 
diation field in the cooling process. In the case of T r « Tk, eq. 
(3) tends to the case considered in Flower et al. (2000). At high 
redshifts, the Cosmic Microwave Background Radiation (CMBR) 
temperature, T r = Tcmbr, may become comparable to the kinetic 
temperature of the gas Tk ; this will significantly affect the HD CF 
at that Tk (see below for more details). 


3 THE COOLING FUNCTION OF HD MOLECULE 

In Fig. 1 we present the variation of Whd with Tk and uh- The 
dotted-line curves correspond to the HD CF, neglecting the colli¬ 
sional ro-vibrational transitions as in Flower et al. (2000). Our re¬ 
sults agree very well with those of these authors. Note that we did 
not consider collisions of HD with He and H 2 because, as it was 
mentioned above, Flower et al. (2000) showed that Whd is essen¬ 
tially a function of Tk and uh only. The solid-line curves in Fig. 1 
are the CFs including ro-vibrational collisional transitions for a 4- 
level vibrational structure of the HD molecule (v = 0,1, 2, 3). One 
sees that the latter CFs depart from the former, as Tk and uh are 
larger. For iih > 10 6 cm~ 3 and Tk « 10 4 K, the calculated value 
of Whd including the vibrational transitions is more than a fac¬ 
tor of fifty larger than the Whd calculated neglecting collisional 
ro-vibrational transitions. Note that in Fig. 1, Whd has been cal¬ 
culated for T r = 2.73 K (the CMBR temperature at present), i.e., 
well within the limit T r « Tk, when the radiation field does not 
participate in the thermal balance of the gas and eq. (3) is reduced 
to the equation for Whd used in Flower et al. (2000). This way 
were able to compare our results with those of these authors. 

The exact HD CF data presented here are avail¬ 
able at http://www.cifus.uson.mx/Personal_Pages/anton/ 
DATA/HD_cooling/HD_cool.html. To facilitate the use of the 
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Table 1 . Polynomial coefficients D; m 



m = 0 

m = 1 

m = 2 

m = 3 

m = 4 

l = 0 

-42.57688 

0.92433 

0.54962 

-0.07676 

0.00275 

l = 1 

21.93385 

0.77952 

-1.06447 

0.11864 

-0.00366 

l = 2 

-10.19097 

-0.54263 

0.62343 

-0.07366 

0.002514 

l = 3 

2.19906 

0.11711 

-0.13768 

0.01759 

-0.000666317 

l = 4 

-0.17334 

-0.00835 

0.0106 

-0.001482 

0.000061926 


HD CF in computational programs we give also a polynomial 
approximation for the HD CF. The approximation is very accurate 
in the ranges of uh and Tk studied here, from 1 to 10 8 cm” 3 , and 
from 10 2 to 2 x 10 4 K, respectively. The approximation is writen 
in the form: 

Loq(Whd) = E DimTlriH D , (4) 

Z ,m=0 

where the coefficients Di m are tabulated in Table 1. 

From Fig. 1 one clearly sees how Whd increases at high den¬ 
sities after including the collisional ro-vibrational level transitions. 
In the case of low densities, the contribution of these transitions 
to the cooling of HD is less important. Nevertheless, at high tem¬ 
peratures Whd is still larger than in the case when these transi¬ 
tions are omitted (pure rotational transitions shown by dashed line 
at Fig. 1). Therefore, the analytical approximation commonly used 
for the low density limit (see e.g., Galli & Palla 2002) should be 
slightly modified. The approximation given in Galli & Palla (2002) 
is based on the calculations for only two low rotational transitions 
(J' J) = (10) and (J'J) = (21), because the collisional probabil¬ 
ities of these transitions are rather large as compared with those of 
higher rotational levels. However, it should be stressed that when 
the vibrational structure is taken into account, then there are other 
ro-vibrational transitions with comparable probabilities of colli¬ 
sional excitation. Based on our results, we suggest the following 
modified approximation for the low-density limit HD CF: 

Log (Whd) = -42.45906 + 21.900837), - 10.1954T 2 + 

2.197887) 3 —0.172867)?. (5) 

This aproximation can be applied for the gas density iih up to 
10 3 - 10 4 cm -3 . 

The effects of the radiation field u in Whd become important 
when T r « Tk (Flower 2000; Flower & Pineau des Forets 2001). In 
this case, the radiative absorption by the HD molecule dominates 
and the second term in eq. (3), as well as the term in the 

probability Wff ,j,, should be considered. The second term in 
eq. (2) can be interpreted as a heating function. At the redshfits at 
which the first baryon objects are expected to form (z ~ 10 — 100), 
the CMB temperature can be similar to the typical kinetic temper¬ 
atures of the small primordial gas clouds. In Fig. 2 we present the 
complete HD CFs (with the collisional ro-vibrational transitions in¬ 
cluded) considering a radiative field with T r = Tcmbr(0)(1 + z) 
for redshifts 2 = 10, 30 and 100 (dotted, dashed and long-dashed 
lines, respectively) and for two densities, uh = 10 and 10 8 cm” 3 . 
One sees that when T r ~Tk, Whd falls dramatically. The second 
term in eq. (3) actually comes to dominate and the HD molecu¬ 
lar lines act as heating sources. This will inhibit the collapse of a 
gas cloud. Therefore the minimum temperature of the primordial 
collapsing clouds is limited by the CMBR temperature (or by the 
heating due to other radiation fields) and not by HD (and H 2 ) cool¬ 
ing. 

We note that the steady state assumption for the HD popu- 
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Figure 1 . Cooling function of the HD molecule, calculated from nn = 1 to 
10 8 cm” 3 in powers of ten (from bottom to top). Radiation temperature T r 
was set equal to 2.73 K (T r « Tk). Solid curves correspond to the cases 
when collisional ro-vibrational transitions (v = 0,1, 2, 3, and J ^ 8) were 
taken into account, while dotted curves are for the case when these transi¬ 
tions were neglected in the calculations. The latter curves can be compared 
with those of Fig. 3 in Flower et al. (2000). Note that the curves in both 
cases saturate for nn > 10 5 — 10®. 


lation levels (see §2) remains even at the high redshifts considered 
here. The typical time scales for population level changes are deter¬ 
mined mainly by the inverse of the transition Einstein coefficients. 
For the HD molecule transitions these coefficientes imply typical 
time scales much less than the time scales related to variations in 
Tcmbr and/or to the gas clouds collapse in the epochs considered 
here (2 < 100). 


4 DISCUSSION 

We have calculated the HD CF, Whd, for a wide range of ki¬ 
netic temperatures and gas densities by taking into account the 
ro-vibrational structure of the molecule for both the radiative and 
collisional transitions. We have found that, when including the 
collisional ro-vibrational transitions, the cooling efficiency of HD 
molecule is higher than previously reported. The main increasing 
of Whd as compared to previous works is at the high-temperature 
side and for high densities (see Fig. 1); for nn > 10 s cm” 3 , 
the differences in Whd reach an order of magnitude and more at 
Tk ~ 10 4 . Our results in support of recent claims (Flower 2000; 
Uehara & Inutsuka 2000; Flower & Pineau des Forets 2001; Naka¬ 
mura & Umemura 2002; Galli & Palla 2002), suggest that the HD 
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Figure 2. Cooling function of the HD molecule (including collisional ro- 
vibrational transitions) calculated for different values of the radiation tem¬ 
perature, T r = 2.73(1 + z)I\, 2 = 0, 10, 30 and 100 (solid, dotted, 
dashed, and long-dashed lines, respectively), and for njj = 10 and 10® 
cm -3 . 


molecule may play an important role in the thermal balance and 
cooling of primordial gas at high densities. 

In Fig. 3 we attempt to compare the relative contributions to 
the CF of the gas by the Fh and HD molecules. Because the ratio of 
abundances of HD to H 2 is significantly smaller than 1, the Whd 
curves should be reduced by a factor equal to this ratio in order to 
compare the contributions to the cooling of gas by both molecules. 
As it was mentioned above, this ratio for the primordial gas, after 
chemical fractionation, is approximately 1CP 2 — 1CP 3 . In Fig. 3, 
we reduce Whd by a factor of [HD\/[H 2 \ = 10~ 2 ' 4 according to 
Galli & Palla (2002) for 2 = 10. As one can see in Fig. 3, the role 
of HD in the cooling of low-density gas is still unimportant with 
respect to the one of H 2 , while for high densities, HD becomes as 
efficient a coolant as H 2 , being even more efficient at low temper¬ 
atures. The H 2 CF has been calculated according to Flower et al. 
( 2000 ). 

The contribution to the gas cooling of HD found here for high 
gas densities and temperatures larger than ~ 3000 K is close to 
the corresponding contribution of H 2 . This will surely affect cal¬ 
culations of galaxy formation in the regime when dense gas cools 
down from temperatures above ~ 3000 K. This regime may arise 
in photoionization-heated gas or in shock-heated dense material 
in dwarf galaxy-size dark matter halos with masses > 10 8 Mg 
(Bromm & Larson 2004). The presence of an efficient coolant, 
besides that H 2 , able to operate at Tk ~ 10000 K and at high 
densities (as is the case of HD) is probably also relevant to under¬ 
stand the formation of the first globular clusters. Bromm & Clarke 
(2002) have shown that, under the action of an efficient coolant at 
Tk < 10000 K, stellar clusters with masses ~ 10 5 — 10 1 may form 
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Figure 3. Cooling functions of the HD (solid lines) and H 2 (dashed lines) 
molecules for two gas densities, n^ = 1 and 10® cm -3 , as indicated 
in the panel. To be compared, the HD CF was reduced by the factor 
[HD]/[H 2 ] = 10 -2 4 . The H 2 CF is from Flower et al. (2000). 

into small dark matter subhalos that later on are tidally destroyed 
during the violent relaxation of dwarf galaxy-sized halos at 2 > 10 
(see also Weil & Pudritz 2001). 

The cooling efficiency of H 2 or HD molecules increases with 
density until the population levels reachs the local thermodynami¬ 
cal equilibrium at a critical density n c , in such a way that at n c the 
probability of collisional de-exitation becomes equal to the sponta¬ 
neous radiative probability. Beyond this density, the cooling is sat¬ 
urated. For molecular hydrogen, the critical density is rather low, 
n c ,H 2 ~ 10 3 — 10 4 cnF 3 , due to the relatively small values of 
the H 2 Einstein coefficients. In the case of the HD molecule, the 
Einstein coefficients are larger as compared to those of H 2 by 2 — 3 
orders of magnitude. One may estimate the ratio of critical densities 
n c ,HD/n c ,H 2 ~ A-v' J' vj ( HD)/A v tjt vj {H 2 ), which is « 10 2 
for radiative vibrational transitions. The value of ti C} hd is indeed 
around two orders of magnitude higher than n c ,ir 2 - From Fig. 1 
(see also Flower et al. 2000), one sees that tic.hd r * j 10 5 - 10 6 
cm -3 . To conclude, we remark the importance of taking into con¬ 
sideration the HD CF calculated here in the simulations and models 
of the first baryonic objects in the universe. In particular, those pro¬ 
cesses related to dense gas that cool down from temperatures above 
~ 3000 K will be affected by our results. 
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